Nuclear Wave Functions for Spin and Pseudospin Partners 
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Using relations between wave functions obtained in the framework of the relativistic mean field 
theory, we investigate the effects of pseudospin and spin symmetry breaking on the single nucleon 
wave functions in spherical nuclei. In our analysis, we apply both relativistic and non-relativistic 
self-consistent models as well as the harmonic oscillator model. 
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I. INTRODUCTION 

Pseudospin symmetry was suggested in the late sixties 
[0, 0| based on the small energy difference between nu- 
clear energy levels with quantum numbers (n,£,j=i+l/2) 
and (n-l,£+2,j=£+3/2). Pseudospin coupling schemes 
were observed and used to explain different aspects of 
nuclear structure 0, 0, IU, El • I n particular, pseudospin 
symmetry was considered in the context of deformation 
and superdeformation [j| , magnetic moment interpre- 
tation H 0, ^| and identical bands 0, How- 
ever, the origin of this symmetry was unknown. It was 
noticed that, in the relativistic mean field (RMF) cal- 
culations, energy differences between pseudospin part- 
ners were small 0], but it was not until 1997 when 
the origin of this phenomenon was understood |l5| in 
terms of near-equility of the absolute values of scalar 
Vs(r) and vector VV(r) potentials 0. This obser- 
vation stimulated a number of works on this subject 
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Although the origin of pseudospin symmetry has been 
explained, it is still an open question why this sym- 
metry is obeyed so well in non-relativistic calculations. 
This question is particularly important for superheavy 
nuclei and nuclei far from stability where shell effects 
are crucial and a subtle change in interaction may effect 
even in changing magic numbers (sec, e.g., discussion in 
Ref. 

Although spin symmetry appears to be broken since 
spin-orbit splittings are very large, it is possible that the 
spin symmetry breaking may be a dynamic symmetry, 
i.e., the energy levels are not degenerate but the eigne- 
functions preserve the symmetry. In this paper we, there- 
fore, compare effects of both pseudospin and spin sym- 



metry breaking for non-relativistic as well as relativistic 
single- nucleon eigenfunctions for spherical nuclei. 

The structure of the paper is the following. In Sec. [H] 
we briefly introduce the main pseudospin relations con- 
sistent with the Dirac equation. Numerical evidences for 
pseudospin dynamic symmetry in nuclei are analyzed in 
Sec. IIIII while the Sec. lIVl is devoted to the spin symme- 
try. Calculations have been performed for a number of 
spherical doubly-magic nuclei. The results presented in 
this paper correspond to 208 Pb, as a representative case. 
Conclusions are presented in Sec. 



II. PSEUDOSPIN SYMMETRY AND THE 
DIRAC HAMILTONIAN 



A. Pseudospin Conditions on the Dirac 
Hamiltonian 



The Dirac Hamiltonian (% = c = 1) 

H = a-p + Vv{r)+fJV s (r)+f3M. (2.1) 

with external scalar Vg(r) and vector VV(r) potentials, 
vanishing space components and non-vanishing time 
component, is invariant under an SU(2) algebra if the 
scalar potential Vg(r) and the vector potential Vy(r) are 
related up to a constant C ps 0, |3(j : 



V s (r) + V v {r) = C ps 
The pseudospin generators 

s 



S = 



s 



(2.2) 



(2.3) 
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form an SU(2) algebra 

[Si, Sj] = ICijk Sk : 



(2.4) 



where s = U p s U p , U p — cr ■ p is the helicity transforma- 
tion 0], s — <T /2, and a{ are the usual Pauli matrices. 

The operators Si commute with the Dirac Hamiltonian 
H ps satisfying conditions l|2.2|l : 

[Si,H ps }=0, (2.5) 
thus generating an SU(2) invariant symmetry of H ps . 



B. Pseudospin Conditions on the Dirac 
Eigenfunctions 

According to the SU(2) invariant symmetry of H ps , 
each eigenstate of the Dirac Hamiltonian H ps with the 
third component of pseudospin A = 5 nas a partner with 
jl = — g and the same energy, i.e., 



H ps ^(r) = S fc ^(r) 
where fi = ±| is the eigenvalue of S z , 
&*£(r)=£S£(r), 



(2.6) 



(2.7) 



and fc are the other quantum numbers. The eigenstates 
in the doublet are connected by the generators S± : 




(2.8) 



The generators (|2.7|l . (|2.8() do not mix upper and lower 
components of the Dirac eigenfunctions 4>^(r). Since 
the spin operates on the lower components only (see Eq. 
1)2.3(1'). one of predictions of this symmetry is that the 
spatial amplitudes of the lower components of the Dirac 
wave functions are identical in shape [Isj . For spheri- 
cal nuclei this means that the lower components of the 
pseudospin doublets have the same radial quantum num- 
ber n = n [27] ] and the same spherical harmonic rank 
I = £ + 1. Therefore, it is natural to label the doublets 
with their pseudospin quantum numbers (n, £, j — £± h). 
The Dirac eigenfunctions then have the form 



(2.9) 



where is the spin function, Y^(6,4>) is the spherical 
harmonic of rank £, and £j = £ ± 1 for j — £ ± | 

As stated above, the radial wave functions of the lower 
components are equal: 



htj=t+h 



(r). 



(2.10) 



On the other hand, the generators for the upper compo- 
nents depend on the momentum as well as the spin so 



they intertwine spin and space. Therefore, in the pseu- 
dospin symmetry limit, the radial wave functions of the 
upper components 9nij=i±±{ r ) satisfy differential rela- 
tions El: 



D.-, 



( r ) 9fUj=l-i{r) = 



(2.11) 



where 



(2.12) 



t+2 



In the non-relativistic limit, g fl ij{r) is associated with 
the single particle wave function while vanishes. 
Relativistic mean held calculations show that indeed 
faii(r) is small compared to D^g^Jr). The factor is 



' ntj \ 

roughly six as we shall see below. 



Dirac Conditions on the Pseudospin Doublet 
States 



Pseudospin symmetry relates upper components in the 
doublets to each other and lower components to each 
other, but pseudospin symmetry does not relate upper 
components to lower components because the pseudospin 
generators <|2.3[) are diagonal. It is, of course, the Dirac 
equation, which relates upper to lower components. 

For spherically symmetric potentials, the Dirac equa- 
tion is reduced to coupled first order differential equa- 
tions in the radial coordinate only leading to |l5j 

D nlj( r )9nlj{ r ) = 

= [2M + V s (r) - V v (r) - E] /^.(r), (2.13) 

Dnij(r)f n ij(r) = 

= \V s (r) + V v (r) + E] 9nlj {r), (2.14) 

where Vs{r) and Vy(r) are spherical potentials and E is 
the binding energy. 

For heavy nuclei, the vector and scalar potentials are 
approximately constant inside the nuclear interior. At 
the nuclear surface the potentials fall rapidly to zero and 
hence outside the nuclear surface both f(r) and g(r) de- 
crease exponentially. Also the nucleon mass is very large 
compared to the binding energy. In the nuclear interior 
Vs(r) — Vv(r) «const, hence 



D fiij( r )9nij( r ) ~ A /r^( r )> 



(2.15) 



where A = 2M + Vg(r) — Vy(r) — E w 6 fin 1 in our 
calculations. Notice that in the pseudospin symmetry 
limit f nij (r) = f Ai (r) and, therefore, Eqs. (l2~m (l2~T4l 
are consistent with Eq. I|2.11|) . 



3 



III. PSEUDOSPIN DYNAMIC SYMMETRIES 

The relation l|2.11|) is strictly fulfilled only under con- 
dition (23, V<? 4- Vv = (7 r - |2S|. Therefore, compar- 
ing the differences between D^-^i (r)g^j_i_i (r) and 

Dnij = £+i( r )9nej = e+i( r )i one can learn about the pseu- 
dospin symmetry breaking effects. The differential rela- 
tions (|2.11|) have been checked previously only for the 
RMF approximation of a relativistic Lagrangian with 
zero range interactions |28| . The pseudospin breaking ef- 
fects have been studied also by taking the integral form of 
Eqs. (|2.11[l but integral relations depend on the boundary 
conditions and hence are less general than the differential 
relations j2(J. 

In this work, we investigate the pseudospin breaking 
effects for spherical double-magic nuclei by carrying out 
three type of calculations. First, we use the standard har- 
monic oscillator (HO) wave functions. Secondly, we per- 
form non-relativistic self-consistent Hartree-Fock (HF) 
calculations with the SLy4 Skryme force [3l| . Finally, we 
perform relativistic mean field (RMF) calculations using 
the Lagrangian with the NL1 parameter set 





-0.1 



-0.1 

^ o 

-0.1 
-0.2 

. 1 



-0.1 
-0.2 



(a) 


1 

i y 








i f 
i / 
i / 
// 
// 
// 
// 
// 


HF - 


X ' / 

\ '/ 
\ '/ 

\ V 

\=/ RMF 


/ ->,. 




\\ 
\\ 

*\. 
o 


/ S~ 
t / 
I / 
I / 
I / 
i / 
i / 
if 
t/ 


HO 




2 


4 6 


8 


2 4 6 8 


■(b) 


/ 

/ *\ 

-1 V 




s. /\>— 






HF 


^ RMF 












3 S « 


HO; 


\^ 6f ; 


2 


4 6 


8 


2 4 6 8 



r(fm) 



A. Comparison within the harmonic oscillator 
model 



For the spherical harmonic oscillator potential (or 
spherical Nilsson potential) we take the analitycal form of 
wave functions with an oscillator frequency hto = 41/A 1 / 3 
[3~i|. Then, one can express D^A^g^Ar) defined by 
Eq. i(2~TTJl as: 



(-l) a (2fi + £-3/2- a)x a 
{h - 1 - a)\a\T(a + 1+3/2)' 



X r )9nh=l+d r ) 



Dnij=e+± V' 1 )y-hij=e+i 



\A (-l) a (2n-2- a)x a 
X ho (n-\-a)\a\T{a + 1 + 3/2)' 

where the envelope function is 



Xn~M) 



2{2v)^(n-iy.i_ x 



and 



\ r(n + £-l/2) 3 



x = r 2v, v = 



x 2 e 



/2 



2h z 

Expressions (|3.1|l can be presented as products 

D ni 3 ( r )9nl.j( r ) = XnZ( X ) P nij( X ) 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



FIG. 1: Numerical check of identities lETTO . fTWl . 

Comparison between D s£j = £_I<7fi£,- = f_I (dashed line) and 

Dfg j= i + ig^i^t+i (solid line) for 1/ (a) and 2d (b) pseu- 
dospin doublets in 208 Pb obtained in different methods: HF, 
HO and RMF. The plot labeled '6f ' shows 6 times scaled lower 
components of the RMF wave function, see Eq. (12.151 and 
related discussion. 



of the common envelope function Xfii( r ) ano - a polynomi- 
als P^ijix) with power expansion coefficients A a (n,£, j): 

n-l 

PnijW =^(-l) a A(n,^,i) x a . (3.5) 

As seen from Eq. i|3.5|) these polynomials are of the same 
order (h — 1) independent of j, whereas the original har- 
monic oscillator eigenfunction with j = £ — ~ involves a 
polynomial of order h while the harmonic oscillator eigen- 
function with j = I + -| involves a polynomial of order 
h — 1 in x. 

As an example, in the lower left corners of Fig. ^ 
(a) and (b) we compare D^ j= ^ + i(r)g-2 j= 2 + i(r) (dashed 

line) with D^ j= i_ i (r)5 fl ^- = ^ _ i (r) (solid line) using ex- 
pression i|3.1|) or its equivalent l|3.4|) . It is seen that 
^■ = l + |(r)5 ft | i= | + i(r) « Z?^^„i(r)(7^ =i _i(r) in 
the whole range of r considered. 

Systematic calculations of many states and nuclei 
have shown that relation 12.11fl holds better as n in- 
creases or I decreases. In order to understand this prop- 
erty we can use the analytical results from Eq. 
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B. Comparison within self— consistent models 

In Figs. □ (a) and (b) we plot D nij=l+i (r)g nij=l+i (r) 
(dashed line) and D-^ = ^_i(r)g n i j= i_i(r) (solid line) us- 
ing the HO model, the non-relativistic HF approxima- 
tion, and the RMF approximation. 

Comparing the non-relativistic and relativistic mean 
field results we see that the agreement is comparable 
or slightly better than the harmonic oscillator results. 
Therefore, the pseudospin symmetry relations are not 
only approximately valid for the relativistic mean field 
eigenf unctions but also for the non-relativistic HF eigen- 
f unctions. Hence we seem to have pseudospin dynamic 
symmetry, that is, the energy levels are not degenerate 
but the eigenfunctions preserve the pseudospin symme- 
try. 

In order to confirm that the radial wave functions of 
the lower components are approximately equal within a 
doublet we also plot them in the case of RMF calcula- 
tions in the lower right corner of Figs.^( a ) an d (b). We 
multiply these wave functions by a factor of 6 in order 
to be comparable to the upper components as suggested 
by Eq. I|2.15[l . Indeed the amplitudes of the lower com- 
ponents are approximately equal |18| . 

IV. SPIN SYMMETRY AND THE DIRAC 
HAMILTONIAN 



FIG. 2: Structure of D-^g-^ for j=l—\ (dashed line, grey 
points), j=(.+ \ (solid line, black points) obtained in the har- 
monic oscillator model for different pairs n, 1 in 208 Pb. The 
first column shows absolute values of coefficients A a (Eqs. 
13.41 13. 5(l ) as a function of a. The second column shows P h i, 
third column - envelope functions Xhf both in terms of di- 
mensionless variable x, defined in equation 13. 3H . 



For this purpose plots of P^ix), A a (n,£,j) are pre- 
sented in Fig. |21 It is seen that P%j[j—i_i.{x) and 
(x) are overlapping more with n increasing. 



While the quantum number n is responsible for the 
similarities between polynomials, the envelope functions 
Xni( r ) strongly depend on I. The higher I is the 
broader is the 'bell' of Xs/( r )- Consequently, when I de- 
creases, the differences between D^._^ + i{r)g^_i + x{r) 

and D~jj_£_ 1 (r)</ fi j._£_ 1 (r) are reduced in the region 
where f^ = £ + i(a:) and P r ^ j= i_i(x) differ. In gen- 
eral, when h 3> £, then A a (n,£,j — £ — |) k 
A a (h,£,j = £+\) and hence D- ij=i _± (r)p fi ^. = |_i(r) « 
D *ij=l+l( r )9iUb=l+x(r). 



A. Spin Conditions on the Dirac Eigenfunctions 

The Dirac Hamiltonian is invariant under an SU(2) 
algebra if the scalar potential, Vs(r), and the vector po- 
tential Vv(fl), are related |15l I3f3l l35| : 



V s (r)-V v (r)=C s 



(4.1) 



where C s is a constant. Hence spin symmetry can occur 
for very relativistic systems like quarks in a meson where 
both Vs(r) and Vy{r) are large 
The spin generators 



S = 



s 
s 



form an SU(2) algebra 

[Si, Sj] — itijk Sfc, 



(4.2) 



(4.3) 



and commute with the Dirac Hamiltonian H s satisfying 
conditions 14.111 



[S u H s } = 0. 



(4.4) 



Thus the operators Si generate an SU(2) invariant 
symmetry of H s . Therefore, each eigenstate of the Dirac 
Hamiltonian H s has a partner with the same energy, 



H s $y r ) = £7**jL(r), 



(4.5) 
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where k are the other quantum numbers and fi = ±5 is 
the eigenvalue of S z , 

S z $^(r) = (i $^(r). (4.6) 

The eigenstates in the spin doublet will be connected by 
the generators S±, 



S± *L(r) 



(4.7) 



The generators <|4.6[1 , (|4.7|) do not mix upper and lower 
components. Since the spin operates on the upper com- 
ponents only, one of predictions of this symmetry is that 
the radial wave functions of the upper components of the 
Dirac eigenfunctions are identical. In the spherical sym- 
metry limit the Dirac eigenfunctions then have the form 



<I> S 



(r) 



9 m(r) [yW(0,0 x ]« 
ifntj{r) [Y^{0,4>) x$ 



(4.8) 



where ij = I ± 1 for j = I ± \. 

The spin generators 14.2|l are related to the pseudospin 

generators by S = 75S75 where 75 = ( ^ ) ' ^ nere ~ 

fore, the conditions are the same as for the pseudospin 
except that now n — > n,£ — » £, g(r) — > / (r), /(r) — > ^(r): 

= D T Uj=(+l M/„^ =£+ 1 (r) (4.9) 

9nij=e-i( r ) = 9n£j=i+±( r )- ( 4 - 10 ) 



and 



For finite nuclei C s — in Eq. 14.1J1 because the po- 
tentials go to zero for large r. Consequently, equal- 
ity Vs(r) — Vy(r) emerges in the spin symmetry limit. 
Therefore, we do not expect spin symmetry to be con- 
served in nuclei since it is known that Vs(r) and Vy(r) 
are both large and of opposite sign. 



B. Spin breaking for different models 

In Figs. 13 (a) and (b) we plot the upper compo- 
nents gnij^t+i (r) (dashed line) and g ntj=l ^i{r) (solid 
line) using the HO model, the non-relativistic HF ap- 
proximation, and the RMF approximation. The Nilsson 
model (HO) shows perfect agreement of course since it 
has a constant spin-orbit potential. However even the 
self-consistent non-relativistic and relativistic mean fields 
show very little difference between eigenstates of the spin 
doublets. 

In the lower right-hand part of Figs. (a) and 
(b), we compare Klj=e-^ D n ej=f-j s ( r ) fnij=i-±( r ) with 
iD ntj=i+ i(r)f nej=e+ i(r) where the factor of 

= (Vs{0) + Vv(0) + Enij)" 1 scales the expression to 
be comparable in magnitude to the upper components, 
according to equation i|2.14|l . The agreement for these 
differential relations is also very good. 
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FIG. 3: If (a) and 2d (b) spin partners' wave functions 
of 208 Pb obtained in HF, HO, and RMF calculations. The 
plot labeled 'f shows the scaled values of D nl - =l _ig nl j =l _i 
(dashed line) and D nt:j=e+ ig nt:j=e+ i (solid line). See text for 
details. 



V. SUMMARY AND CONCLUSIONS 

In the pseudospin symmetry limit the radial wave func- 
tions of the upper components of pseudospin doublets 
satisfy certain differential relations. We demonstrated 
that these relations are not only approximately valid 
for the relativistic mean field eigenfunctions but also for 
the non-relativistic Hartree-Fock and harmonic oscillator 
eigenfunctions. Generally, we expect them to be approx- 
imately valid for eigenfunctions of any non-relativistic 
phcnomcnological nuclear potential that fits the spin- 
orbit splittings of nuclei. Likewise in the spin symmetry 
limit the radial amplitudes of the upper components of 
the Dirac eigenfunctions of spin doublets are predicted 
to be equal and this is approximately valid for both non- 
relativistic and relativistic mean field models. Also the 
spatial amplitudes of the lower components of the Dirac 
eigenfunctions of spin doublets satisfy differential rela- 
tions in spin symmetry limit and these relations are ap- 
proximately valid in the relativistic mean field model. 

Hence we seem to have both spin and pseudospin dy- 
namic symmetry; that is, the energy levels are not de- 
generate but the eigenfunctions well preserve both sym- 
metries. For both of these symmetries to be conserved 
both the vector and scalar potentials must be constant. 
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Of course this is not true. However, for heavy nuclei 
this is approximately true in the nuclear interior and ex- 
terior. Only on the surface are the potentials changing 
rapidly. This leads to a dynamic symmetry for both spin 
and pseudospin. The spin-orbit splittings are determined 
by d[Vv(r) — Vs(r)]/dr while the pseudospin-orbit split- 
tings are determined by d[Vv(r) + Vs(r)]/dr. Therefore, 
the energy splittings for spin doublets are larger than for 
pseudospin doublets because Vy (r) — Vs (r) changes more 
rapidly on the nuclear surface than Vs(r)+Vv{r) because 
V v {r)-V s (r) > \ V s (r) + V v (r) \ in the interior and both 
go to zero in the nuclear exterior. 
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